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Abstract. These proceedings show the preliminary results of the dijet cross sections and 
the dijet longitudinal double spin asymmetries All in polarized proton-proton collisions at 
y/s = 200 GeV at the mid-rapidity \rj\ < 0.8. The integrated luminosity of 5.39 pb _1 collected 
during RHIC Run-6 was used in the measurements. The preliminary results are presented as 
functions of the dijet invariant mass My y The dijet cross sections are in agreement with next- 
to- leading-order pQCD predictions. The All is compared with theoretical predictions based on 
various parameterizations of polarized parton distributions of the proton. Projected precision 
of data analyzed to date from Run-9 are shown. 



1. Introduction 

The jet production rate in polarized proton-proton collisions is sensitive to the polarized gluon 
distribution Ag(x,Q 2 ) of the proton. The polarized gluon distribution is of particular interest 
in the proton spin physics because the first moment of this distribution is the fraction of the 
proton spin carried by the gluon spin, AG. 

The AG is one of the four terms in the decomposition of the proton spin in the infinite 
momentum frame: 

± = ^ + AG + L q + L g , 

where S, L q , and L g are the contributions from quark spin, quark orbital motion, and gluon 
orbital motion. The quark spin contribution E has been measured using polarized deep inelastic 
scattering (pDIS) data combined with neutron and hyperon f3 decay data [1]. 

One of the primary objects of the spin physics program at RHIC (RHIC-Spin) is to determine 
AG by using polarized proton-proton collisions. An advantage of using proton-proton collisions 
is that gluons participate in high-pr events at the leading order. However, it is challenging to 
determine the kinematics of parton-level interactions, which are desirable quantities to determine 
in an experimental study of the structure of the proton. 



In order to determine the kinematics of leading-order parton-level interactions in proton- 
proton collisions, the momenta of both outgoing partons are needed: 



Xl = ^(e+^ 3 + e +2/4 ), x 2 = ^(e~ y3 + e^ 4 ), 

V s V s 

where the subscripts 1, 2 indicate the incoming partons of the hard interactions, 3, 4 indicate 
the outgoing partons. The momenta of outgoing partons can be estimated by observing two 
final state objects in events such as dijets and photon-jets. Furthermore, the invariant mass and 
average pseudorapidity of the two objects are sensitive to the kinematics according to: 

M 33 V X l X< 2 S i o = hi . 

2 2 X2 

These proceedings show the preliminary results of the longitudinal double spin asymmetries 
All of the dijet production as a function of the dijet invariant mass My: 

All - an + aU - 

The arrows (ti) indicate the orientations of the helicities of the colliding protons. All is sensitive 
to the polarized parton distributions. In fact, in the framework of QCD factorization, All can 
be written as follows: 



All = — 



J2j dx iJ dx 2 Af i (x 1 ,Q 2 )Af j (x 2 ,Q 2 )aLLH^se") 

j dXl j dx 2fi( X liQ 2 )fj( X 2iQ 2 )°( C080 *) 



where Afi is the polarized distribution of the parton i, fi is the unpolarized one, a is the parton- 
level cross section, clll is the parton-level longitudinally double spin asymmetries, and i and j 
run over quark flavors and gluons. 

The proceedings also show the preliminary results of the dijet cross sections as a function 
of the dijet invariant mass My. The cross sections are measured to support the theoretical 
framework which relates polarized parton distributions and measured All by showing the 
unpolarized equivalent of the same framework can quantitatively relate the unpolarized parton 
distributions and measured cross sections. 



2. STAR Detector 

STAR, the Solenoidal Tracker At RHIC, was built to measure wide varieties of nuclear 
interactions in high energy heavy ion collisions and polarized proton collisions [2]. Figure 1 
shows the cross sectional view of the STAR detector. The detector subsystems particularly 
relevant to the measurements presented in these proceedings are the Time Projection Chamber 
(TPC), the Barrel Electromagnetic Calorimeter (BEMC), and the Beam-Beam Counters (BBC). 

The Time Projection Chamber (TPC) [3] is the primary tracking system of the STAR 
detector. It has a cylindrical shape operated within a solenoidal magnetic field of 0.5 T and 
provides the momentum measurements over a range from 100 MeV to 30 GeV. Its acceptance 
is \rj\ < 1.8 with full azimuth. 

The Barrel Electromagnetic Calorimeter (BEMC) is the primary calorimeter at mid-rapidity 
[4]. It is a cylindrical annulus which surrounds the TPC. The BEMC covers \rj\ < 1 with full 
azimuth and has a depth of about twenty radiation lengths (20Ao) at 77 = 0. The BEMC has 
4,800 towers in total and each tower covers A77 x Atp = 0.05 x 0.05. 




Figure 1. The STAR detector 

The Beam-Beam Counters (BBC) [5] are scintillator annuli of hexagonal tiles. Their 
acceptance is approximately 3.3 < \rj\ < 5.0. The BBCs are used to trigger minimum bias 
(MINB) events. The MINB condition is a coincidence between the east BBC and the west BBC. 
The cross section of the MINB events is <tmb = 26.1 ± 2.0 mb [6]. The MC simulation estimates 
that 87% ± 8% of non-singly diffractive collisions result in a MINB trigger [5] . 

3. Event Selection 

The BJP1 (Barrel Jet Patch 1) trigger was primarily used in the measurements. This 
trigger requires a minimum transverse energy E^ deposit in a patch of calorimeter towers 
(A77 x Acp = 1 x 1) as well as the MINB condition. The threshold required for offline 
analysis was 10.8 GeV, which was above the trigger turn-on, to ensure high trigger efficiency. 

The BBC coincidence has some allowed time difference. This difference is measured as 4-bit 
values called timebin, which roughly corresponds to the vertex position of the events. To select 
events close to the interaction point (IP), events are required to be in a specific timebin. The 
vertex distribution of the events in this timebin has approximately a Gaussian distribution with 
the mean -19 cm and the standard deviation 30 cm. About 26% of the MINB events were in 
this timebin. In addition, events are required to have a reconstructed vertex. 

4. Jet and Dijet Definition 

Jets can be defined at three different levels: the parton level, the hadron level, and the detector 
level. In an experiment, jets are reconstructed at the detector level, whereas perturbative QCD 
calculations predict jet productions at the parton level. MC simulated events are used to evaluate 
the effects of the transitions between different jet levels as jets can be reconstructed at all three 
levels in MC simulation. 

Jets are defined as collections of four-momenta selected by the mid-point cone jet-finding 
algorithm [7] with the cone radius 0.7 and split/merge fraction 0.5. Four- moment a of jets 
are the four-vector sum of the four-momenta that define the jets. Four-momenta that compose 
detector- level jets are constructed from charged tracks in the TPC and energy deposits in BEMC 
towers. Tracks are assumed to have the mass of a charged pion (139.75 MeV), and towers are 
assumed to be massless. In order to avoid measuring the same charged particles twice both in 



the TPC and in the BEMC, if a track points to a BEMC tower, the energy that a MIP would 
leave in the tower is subtracted from the tower energy. 

In order to reject the beam- gas background, the neutral energy ratio Rt, the fraction of 
jet transverse energy E^ reconstructed from energy deposits in the BEMC, are required to 
be smaller than particular values that depend on jet p^\ Rt < 1.0 (5 < pr < 17.31 GeV), 
R T < 0.99 (17.31 < p T < 21.3 GeV), R T < 0.97 (21.3 < p T < 26.19 GeV), and R T < 
0.90 (26.19 GeV < p T ). 

Dijets are defined as the two leading-px jets of events. Dijets are required to contain at least 
one trigger jet, a jet which caused the BJP1 trigger. If only one jet is a trigger jet, the jet is 
called the same side jet. If both jets are trigger jets, the jet that triggered with higher E^ is the 
same side jet. The other jet is called the away side jet. Dijets are required to have balanced p^: 
0.73 < p^ way /p^ me < 1.1 because px-balanced dijets are more likely to carry momentum closer 
to that of parton level than do unbalanced jets. 

Asymmetric p^ cuts (max(^x) > 10.0 GeV and min(^x) > 7.0 GeV) are used because a NLO 
pQCD calculation has little prediction power of dijets cross sections and All with symmetric 
Pt cuts. The —0.8 < 77 < 0.8 cuts are applied for the detector acceptance. The 1 773 — 774 1 < 1.0 
cut is necessary for both jets to be in the acceptance at the same time. The A<p > 2.0 cut is 
applied to select back-to-back dijet events. 

5. MC Simulation 
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Figure 2. The data-MC comparison of the dijet kinematic distributions, (left) The invariant mass 
Mjj distributions, (top right) The average pseudo-rapidity r] distributions, (bottom right) The pseudo- 
rapidity difference A77 distributions. The MC yields are scaled so that the yield becomes the same as the 
data. 



The events are generated by the Pythia 6.410 event generator [8] with the CTEQ5L parton 



distributions [9] using parton between 3 and 65 GeV. The detector responses to the events 
are simulated with a GEANT3 [10] based STAR detector simulation program. 

In the MC simulation, jets are reconstructed at all three jet levels using the same jet finder 
as in the data. The detector-level jets are defined in the same way as in the data. The hadron- 
level jets are collections of final state particles in the event generator, while the parton-level jets 
are composed of outgoing partons of the hard interactions and the radiation from the outgoing 
partons. 

The MC events well reproduce the data. Fig. 2 shows the data-MC comparison of the dijet 
kinematic distributions: the invariant mass My, the average pseudo-rapidity 77 = (773 + t?4)/2, 
and the pseudo-rapidity difference A77 = 773—774 distributions. 

6. Cross Section Measurement 

The dijet cross sections are estimated for each My bin with the formula: 

d 3 a _ 1 1 1 

dMjjd7/ 3 d7/4 ~ JCdt AMjjA7/ 3 A7/4 C ' 

J is the dijet yields at the detector level. C is the correction factors which correct the dijet yields 
to the hadron level. The correction factors C are estimated from the MC events as bin-by-bin 
ratios of the dijet yields at the detector level and at the hadron level. AM^A^A^ normalizes 
the cross sections to per unit space volume. J Cdt = 5.39±0.41 pb _1 is the integrated luminosity 
measured using the BBCs. 

The major systematic uncertainty is due to the uncertainty on the jet energy scale (JES). 
Because the My dependence is steeply decreasing, the uncertainty of the cross section is very 
sensitive to systematic uncertainty on the JES. 

The track portion of the jet energy has 5.6 % of systematic uncertainty. To evaluate the 
effect of this uncertainty on the cross sections, the cross section was reevaluated with the 5.6 % 
variation of the track portion of the jet energy. 

Energy deposits in the BEMC towers have 4.8 % of systematic uncertainty. The effect of this 
uncertainty was evaluated by varying the BEMC tower energies. After the energies were varied, 
the offline trigger thresholds were reapplied and the jet finding algorithm was rerun. 

Systematic uncertainty due to the correction for the pile-up and the timebin selection are 
estimated to be small compared to the JES uncertainty. The cross section has 7.6% of correlated 
systematic uncertainty due to the uncertainty of the integrated luminosity. 

The dijet cross sections are calculated by next-to-leading order perturbative QCD with the 
CTEQ6M [11] parton distributions. Jets are defined by a cone jetflnding algorithm with cone 
radius 0.7. Both the renormalization scale and the factorization scale are \i = My. The scale 
uncertainty is calculated by varying the scale from 0.5My to 2My. 

The effects of the hadronization and the underlying events were evaluated by using the MC 
events. The correction factors Chad were obtained for each My bin as the ratios of the cross 
section at the hadron level and at the parton level. The systematic uncertainty on Chad was 
calculated by varying the cone radius from 0.6 to 0.8. 

Figure 3 shows the preliminary results of the dijet cross section measurements. The measured 
dijet cross sections are well described by the theory. This indicates that measured All as 
well can be interpreted in the same theory and suggests ways to constrain the polarized gluon 
distributions from dijet All- 




Figure 3. (left) The dijet cross sections compared to theoretical predictions. The systematic uncertainty 
does not include 7.68% of uncertainly due to the uncertainly on the integrated luminosity, (right) The 
ratios: (data - theory) /theory. The theory includes the NLO pQCD predictions and the corrections for 
the effects of hadronization and underlying events. The ratios are taken bin by bin. 



7. Longitudinal Double Spin Asymmetry All Measurement 

The dijet longitudinally double spin asymmetries All were measured as the ratios of the spin 
sorted dijet yields with the corrections for the relative luminosity and polarizations: 

A _ E PyPB{(N n + N U ) - R(N n + N lt )} 

LL E P$P£{(N n + n u ) + R(n u + N lt )Y 

The arrows (ti) indicate the orientations of the helicities of the proton beams. The relative 
luminosity R — (C^ + C^) / + was measured using the BBCs. The relative luminosity 
varied between approximately 0.9 and 1.1. The polarizations (PyPs) was measured by the pC 
CNI polarimeter and the polarized H jet polarimeter [12]. The average polarization for the 
Yellow beam and the Blue beam of RHIC were Py = 59% and Pb = 56%, respectively. The 
figure- of -merit = PyP^C for this measurement is 0.59 pb" 1 . 

Four false asymmetries, which should vanish, are measured for a systematic check of the data. 
Two single spin asymmetries and two wrong-sign spin asymmetries were consistent with zero 
within the statistical uncertainties. 

The largest systematic uncertainty is the trigger bias, which is the uncertainty in the changes 
of All from the parton level to the detector level. This was evaluated using the MC events with 
several polarized parton distributions which are compatible with the current experimental data: 
DSSV, GRSV std, and the GRSV series with AG from -0.45 to 0.3 [13], [14]. In this evaluation, 
in order to calculate All with the unpolarized event generator Pythia, the MC events were 
weighted by the products of the spin asymmetries of the parton distributions and parton-level 
cross sections: (A/i(^ H , Q? )Af 2 (x 2i , Q\ )/fi(x lh Q? )/ 2 (^2i 5 Q\ )) ■ a LL (cos#*). 



Figure 4 shows the preliminary results of dijet All measurements. It can be seen that 
the results are consistent with the next-to-leading perturbative QCD prediction of the DSSV 
polarized parton distributions. The results are also consistent with GRSV zero scenario, and 
larger values of Ag in GRSV are disfavored. 
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Figure 4. The double longitudinal spin asymmetry All for the dijet production as a function of dijet 
mass Mjj. The vertical bars on the data points indicate the size of the statistical errors. The horizontal 
bars on the data points indicate the bin widths. The data points are plotted at the mean values of Mjj 
of the events in the bins. Predictions of next-to-leading perturbative QCD with various models of the 
polarized gluon distributions are shown. 



8. Projected Sensitivity 

Data collected during 2009 at RHIC represents a considerable increase in sensitivity to All for 
dijet production. Approximately 22 pb" 1 were recorded with an average polarization of 59% in 
both beams. Data corresponding to a figure-of-merit of 0.96 pb" 1 have been processed using 
the TPC and BEMC to measure dijets at mid-rapidity. With the additional statistics, the data 
can be divided into different pseudorapidity acceptances, which provides sensitivity to the ratio 
of xi/x2. The invariant mass distribution for different pseudorapidity acceptances is therefore 
sensitive to different x\, X2 phase space and allows extraction of constraints on the shape of 
Ag{x). 

The two interesting divisions are when the two jets have the same sign pseudorapidity and 
when they have opposite signs. The statistical precision of the data analyzed to date can be 
seen in Fig. 5, which is improved over comparable figure-of-merits from previous years due 
to improvements in trigger efficiency. An expected increase in the figure-of-merit by a factor 
of between 1.5 and 2.0 will provide further improvement of the sensitivity. Various studies to 
understand systematic uncertainties are underway. 

9. Summary 

These proceedings showed the preliminary results of the dijet cross sections and dijet All in 
polarized proton-proton collisions at — 200 GeV from Run-6. These results agree well 
with pQCD calculations and provide constraints on the gluon polarization Ag(x) in the proton. 




Figure 5. The statistical precision of 2009 STAR dijet data analyzed to date. The data has been divided 
into two pseudorapidity regions, which provide different sensitivities to the Bjorken-x phase space. The 
data in the left (right) panel are from when the two jets have the same (opposite) sign pseudorapidity. 



Projected precision from Run-9 data in multiple acceptances shows that STAR will be able to 
add new constraints to the shape of Ag(x). 
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